Quantum Fluctuations Affect the Critical Properties of Noble Gases 
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Interacting argon atoms are simulated with a recently developed quantum Langevin transport 
treatment that takes approximate account of the quantum fluctuations inherent in microscopic 
many-body descriptions based on wave packets. The mass distribution of the atomic clusters is 
affected significantly near the critical temperature and thus it may be important to take account of 
quantum fluctuations in molecular-dynamics simulations of cluster formation processes. 
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Molecular dynamics presents a powerful tool for elu- 
cidating the features of mesoscopic systems The 
present communication draws attention to the possible 
importance of quantum fluctuations in such treatments. 
To illustrate this issue, we focus on the cluster mass dis- 
tribution for argon atoms in thermal equilibrium. 

Entities interacting via van der Waals forces, such as 
atoms of noble gases, form clusters exhibiting enhance- 
ments at magic mass numbers associated with especially 
compact geometric configurations Several studies 

based on molecular dynamics have been devoted to the 
exploration of the mass spectrum of argon clusters in 
various scenarios (see for example [§-0). However, the 
detailed formation process of the clusters, and their ob- 
served relative abundance, is not yet fully understood 
and may well depend on the specific experimental con- 
ditions. Because of this delicacy, we examine here the 
possible role of quantum fluctuations. Such fluctuations 
are in principle present in the the classical description be- 
cause it is based on the assumption that the many-body 
wave function can be regarded as a product of single- 
particle wave packets, which are not energy eigenstates. 

Specifically, our discussion invokes a recently devel- 
oped treatment in which this inherent effect is emulated 
by a quantal Langevin force augmenting the usual classi- 
cal molecular-dynamics equations of motion [^|J|]. That 
treatment was developed in the context of nuclear frag- 
mentation processes and we first briefly describe its adap- 
tation to atomic systems. Subsequently we use the model 
to calculate the cluster mass spectrum for thermal ensem- 
bles of argon atoms and discuss how the results compare 
with the usual classical simulation. 

The interaction between two atoms i and j is described 
by the usual Lennard-Jones potential. 
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where the values a = 3.405 rA and e = 119.8 K are ap- 
propriate for argon atoms and they provide convenient 
scales for distance and energy, respectively. (For numeri- 



cal convenience, the potential is truncated at ri 
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in the potential occurs at an interatomic separation of 
1.1225 a. Furthermore, the strength of 
the interaction is characterized by the parameter e which 
equals the magnitude of the potential at its minimum, 
Vo = —e, and provides a convenient energy unit. 

The simulation of a classical fluid in thermal equilib- 
rium can be accomplished by augmenting the molecular- 
dynamics equations of motion by a stochastic term rep- 
resenting the coupling to a thermal reservoir. 
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Here fi^— J2j the force exerted on the atom i 

by all the other atoms j ^ i via the employed Lennard- 
Jones potential (Q). Furthermore, the Langevin force 
C, emulating the coupling to the reservoir is a standard 
white noise with the correlation function 
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where / is the 3x3 unit matrix. The rate v deter- 
mines the strength of the coupling and should be cho- 
sen to be sufficiently small that its value does not af- 
fect the results. Finally, the matrix M has the elements 
Mij = {5ij — 1/N)I, and so the Langevin force does not 
affect the overall momentum P = J^i Pi- 

When the equations of motion are propagated by a 
simple leap-frog method, the increment in the atomic 
momenta can be calculated as Pj Pi + Ap^ with 



Ap, - f^At - vAtp^ + ViJAtfx, 
^{vAtY,P,~^/^Y.^,) 
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which has no effect on the results.) The minimum 



where Xi s-re a random vectors whose components have a 
normal distribution, -< Xi >~— and -< XiXj >~— '^SijI. 

In the quantal treatment, each individual atom is de- 
scribed by a Gaussian wave packet with a spatial variance 
equal to Ar^ = h/2mu!, where m is the atomic mass 
and the frequency lo is determined by from the relation 
iv^ = Vq /m. The corresponding "level spacing" is 
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The inclusion of the quantum Langevin force effectively 
modifies the drift term for the intrinsic cluster motion by 
a factor given by PjlO|] 



I(l-e-'-)<l 
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In terms of this reduction factor a, the quantal Langevin 
equation can be written as follows. 



p = / - lyM-" ■ [aip - mV) + mV] + VvTM • C , (8) 

where V is the velocity of the cluster to which the atom 
belongs. The only complication arises from the require- 
ment that the center-of-mass motion of each individual 
cluster retain a classical character. Thus, in order to de- 
termine V, it is necessary to perform a cluster analysis of 
the iV-body system at each time step. This task can be 
easily accomplished by considering two atoms as belong- 
ing to the same cluster if (and only if) their separation 
is less than a certain critical value, taken to be equal 
to 2a. When the drift term is modified correspondingly, 
the expression (^) for the momentum increment is being 
replaced by 



where V„ is the velocity of the cluster n containing i. 

The occurrence of the quantum Langevin force is a 
consequence of the fact that the wave packets are not 
energy eigenstates. That same feature causes another 
complication, namely the need for taking account of the 
thermal distortion of the intrinsic structure of each wave 
packet (the relative weight of each energy component of 
the wave packet is temperature dependent) . This feature 
can be taken into account by subjecting the dynamical 
state of the system to a cooling process before making 
the observation |^,|0|. In the present case, this amounts 
to integrating the equation 
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from T = to T = h/2T, thereby distorting the cur- 
rent dynamical values {r{t),p{t)) into {r' (t) , p' (t)) . As 
a result of this cooling process, the atoms change their 
positions and so the cluster structure may be modified. 
While this is a significant effect in the case of nuclear 
fragmentation |10|, the effect is unimportant for atomic 
systems, since the width of the wave packet Ar is quite 
small (less than 5% of (T in the case of argon). 
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FIG. 1. The mass distribution of argon clusters calculated 
either wtth (solid circle) or without (open triangle) the quan- 
tum Langevin force emulating the quantum fluctuations in- 
herent in wave-packet dynamics. The temperatures are indi- 
cated in units of e and the atomic density is p = 0.025 a~^. 

In order to illustrate the effect of the energy fluctu- 
ations caused by the use of wave packets, we show in 
Fig. |l] the cluster mass distributions obtained for a sys- 
tem of argon atoms held at a specified temperature T. 
These calculations consider 100 atoms confined within 
a cubic box and subjected to periodic boundary con- 
ditions. The temperatures considered range from well 
above to below critical. At high temperatures the clus- 
ter mass distribution falls off rapidly, whereas there is 
a preference for condensation into large clusters at low 
temperatures. When the classical Langevin equation (||) 
is employed, the critical temperature associated with this 
gas-liquid phase transition is approximately T'^' w 0.6e, 
as can be seen from the figure. The inclusion of the quan- 
tum Langevin force, Eq. (||), produces a steeper slope in 
the mass distribution at high temperatures and the phase 
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transition occurs at a lower temperature, T'^" « 0.5e, as 
might then be expected. 

The change of the critical temperature can be under- 
stood quantitatively by noting that it is possible to ex- 
tract an effective classical temperature from the quantal 
Langevin equation (H) . The effective temperature can be 
estimated by means of the Einstein relation as the square 
of the diffusion coefficient divided by the drift coefficient, 

r,ff = ^ = ^ = i?/(i-c-^/^)>r. (11) 

aiy a 

This estimate ignores the correction terms arising from 
the cluster center-of-mass motion, as is justified when we 
are interested in the cluster mass distribution, because 
this observable depends mainly on the intrinsic degrees 
of freedom of the cluster. In addition, the cluster mass 
distribution is not affected much by the cooling process, 
as already mentioned. Thus, it is expected that a calcu- 
lation with the classical Langevin model at the effective 
temperature, T^s, will yield results that are similar to 
those obtained with the quantal Langevin model at the 
real temperature, T. 
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FIG. 2. The cluster mass distribution obtained with ei- 
ther the quantal Langevin model at T = 0.5e (solid circles) 
or the classical Langevin model at the corresponding effec- 
tive temperature T^s = 0.62e (open triangles), at the density 
p = 0.025 cr-^ 

This expectation is indeed borne out, as shown in Fig. 
H where we compare the cluster mass distribution ob- 
tained with the quantal model at T = 0.5e to the result of 
the classical treatment carried out at the corresponding 
effective temperature Tcs — 0.62e. The quantitative sim- 
ilarity between the two distributions is remarkable and 
supports the above discussion. 

The temperature shift depends the system under con- 
sideration and the specific observables extracted. For 
example, the temperature shift is in the opposite direc- 
tion for particles in a harmonic oscillator potential [p|JTl| 
and for nuclear fragmentation jl^. This is because the 
cooling process affects the observables as well as the mod- 



ification of the drift term for these systems. In order to 
illustrate this feature, we consider here the evolution of 
the distorted momentum (i.e. the solution to Eq. (|lC 
which is given by 
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Thus, in the rest frame of the cluster, the distorted mo- 
menta of the constituent atoms are governed by a modi- 
fied Langevin equation. 



p' = C-^/2T J _ ^^]^2 . _^ 



°^^^V^M ■ C ■ (13) 



Since this equation of motion has a classical form (c./. 
Eq. (H)), we can again invoke the Einstein relation and 
extract an effective temperature for the intrinsic cluster 
motion. 



= D/(e^/^- 1) < T . (14) 



It has been shown that calculations with classical molec- 
ular dynamics at this equivalent temperature T'^^ yields 
results that are very similar to the exact quantal results 
for the real temperature T, for non-interacting particles 
in a harmonic potential |^,^ It is to be expected 
that this feature holds as well for the distribution of in- 
trinsic excitation energies of atomic clusters. Therefore, 
it would be interesting to measure the atomic cluster 
mass distribution simultaneously with the intrinsic exci- 
tation energy distribution, since the associated tempera- 
ture shifts are predicted to be of opposite sign (and thus 
they could not be mocked up by a simple readjustment 
of the classical temperature). 

The present treatment takes account of the fact that 
the intrinsic atomic motion in a cluster and its center-of- 
mass motion are of different character. While the latter 
can be considered as classical, the quantal features of the 
former may be significant and lead to modifications of the 
Einstein relation, as we have discussed in the framework 
of wave-packet dynamics. This important point was al- 
ready noted by Ikeshoji et al. in their recent study of 
magic-cluster formation and it might thus be interesting 
to apply the present treatment to the cluster formation 
process in an expanding and cooling noble gas. 

In the present paper we have adapted a recently devel- 
oped quantal Langevin treatment to a system of argon 
atoms in thermal equilibrium. The method was de- 
veloped in the context of nuclear dynamics and takes 
approximate account of the energy fiuctuations that are 
necessarily present when wave packets are used to de- 
scribe the system. The presence of these quantum fluc- 
tuations changes the character of the specific heat from 
classical to quantal and their inclusion by the developed 
method leads to a significant improvement of the statisti- 
cal properties in a number of simple test cases that can be 
subjected to exact analysis 11 1 as well as for finite nu- 
clei [plpl . When incorporated into microscopic dynamical 
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simulations of nuclear collisions, it leads to a significant 
improvement of the calculated fragment mass distribu- 
tion For atomic systems, the effect of quantum 
fluctuations is generally expected to be small, due to the 
relatively large mass of the "particles" . Nevertheless, as 
we have shown above, the inherent quantum fluctuations 
may also play a role in atomic physics. Certainly, as our 
study brings out, it appears that the critical properties of 
noble gases are affected significantly. The quantum fluc- 
tuations may therefore also affect the formation process 
and might affect the outcome of dynamical simulations 
aimed at understanding the observed mass distributions. 
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